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МАТЕМАТИКА 





УДК 519.642.3 


В.М. ДРАГИЛЕВ 


О НЕВЯЗКЕ ИНТЕГРАЛЬНОГО УРАВНЕНИЯ ПЕРВОГО РОДА 
В МЕТОДЕ ПРОЕКЦИИ 


В предшествующих работах для интегрального уравнения Фредгольма первого 
рода с вырожденным ядром развивалась обобщенная (не регуляризующая) вер- 
сия метода проекций, основанная на произвольном выборе базиса в сочетании с 
некими априорными оценками для погрешности решения. Эти оценки зависят от 
невязки, которой ранее пренебрегали. Результаты настоящей работы позволяют 
учитывать невязку и показывают, что, вообще говоря, ее учет является суще- 
ственным. 

Ключевые слова: интегральное уравнение Фредгольма первого рода, метод 
проекции, априорные оценки погрешности. 


Введение. Рассмотрим интегральное уравнение Фредгольма первого рода 
с вырожденным ядром 


(1) 
где и(х)=и(х)+ди(х); ис) =[А4]©); 9(5)е Г.[а,Б] - искомая 


функция (оригинал); ди(х) - погрешность исходных данных; 


Ф„(х)и„(5) - гладкие комплекснозначные функции; системы 


ии) и {Фи линейно независимы на отрезках 


[аб] и [с,а] соответственно. 
Задача (1) некорректна и может решаться методом Тихонова [1]. В 
другом подходе, известном как метод проекций [2], в пространстве 


[0] 


и_/ И Обобщенное 


Г.[а,Б] вводится ортонормированный базис { }, (5) } 


решение строится в виде 
№ 
4, (5) =>, 7,(5), (М<М) (2) 
В= 


при этом коэффициенты о отыскиваются из условия минимизации не- 


вязки в пространстве [.,[с, 4]. В канонической версии метода проекций 


базис образуют собственные функции (СФ) оператора А”А, что делает 
метод регуляризующим [1, 2] и в неком смысле минимизирует погрешность 
решения [3]. На практике более удобной и достаточно эффективной [4] 
может оказаться обобщенная версия метода, в которой базис выбирается 
волевым образом. 
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Введем обозначения: Ну - линейная оболочка функций 


(ор; Р, - проектор из Г[а,6] в Ну; Ау=АР,у; 
4к($) = Р/4($); Их =|9-а»|14|; А = [би 


определяемых на отрезке [а,Б] или [с,@], символы | . | здесь и ниже 








; для функций, 








означают норму в пространстве [.,[а,6] или Г[.,[с,4] соответственно 


(подобное же соглашение примем для символов скалярного произведе- 
ния). Преимущества метода проекций связаны с возможностью оценить 


погрешность решения Йх — ах —а|/“ при априорно известных по- 
грешностях 7у и А. А именно, справедлива оценка [3] 
Як ЗПу +1 +ИмИм, где м = |@ — а» |Л|а|, 7» = |4 - а» |“ 


Ч, (5) - решение (2) задачи (1) с ди = 0 в правой части. Погрешность 




















Й 


Я удовлетворяет оценкам [3]: 


2 
Пы ЗП < 1+ ХыйЙ, (3) 
где У»м- так называемая норма [3], зависящая от выбранного базиса; в 


К [С 
каноническом методе проекции Ум =0. Для погрешности 7] 
имеется оценка [5, 6] 


Пм <УкСмА, (4) 


; им = Аам; С» - число обусловленности М - мер- 








пе 
ной матрицы БВ» эрмитова оператора А’А/Н,—Ну в 


м 
базисе {],(5)}„_!. Множитель У, рассматривается здесь как функцио- 


нал от 4(5). 
Для применения оценки (4) требуется при произвольно заданной 


величине 7], оценить сверху тах У , что и является целью настоящей 
9(5) 


работы. 
Основной интерес представляет случай 7] <<1, так как в про- 


тивном случае восстановление оригинала с малой погрешностью у <<1 
а рйой исключается. При некоторых предположениях из оценки (3) и из 
непрерывности оператора А следует, что при Пу —>0 обе функции 


и(хуим(х) близки к [Аду |(Х), и в таком пределе У, =1. Поэтому в 
первом приближении фактором У, в (4) пренебрегают [4], а оценивание 


' 
этого функционала интерпретируется как учет невязки | = Им | . 


176 


Вестник ДГТУ, 2006. Т.6. №3(30) 15ВМ 5-7890-0368-0 








Основные соотношения. Будем предполагать, что матрица Бу 
наибольшей размерности № = М не вырождена; тогда в силу так назы- 
ваемого принципа промежуточности [7] матрицы БВ» размерностей 


№ <М тоже не вырождены. В этих условиях 7/=Ги Ух>1 при 
№ < М, причем нижняя грань У,=1 достигается для всех ненулевых 


оригиналов (5) е Нм [6]. Представим У, в виде 
Ум =7м7м,  (М<М» (5) 
Бим = ммм ик о =[Аам 9. 


М 
Пусть (” (5)}„_, есть ортонормированный набор СФ оператора 


(1) 


пе = [| 














* М2 м 

АхАу, отвечающих его собственным значениям {(@, )}_1; 
М№ 

гдеб’, - так называемые сингулярные числа (СЧ), упорядоченные по убы- 


М М № 
ванию: О 20. >...0 у >0. Обычным путем (см., например, [5]) 
устанавливается так называемое сингулярное разложение: 


№ 
[4,9] =У`ои (ЧФ (©, (6) 

ПЕ 
где Ф\ (= (9^) Ао) - функции, ортонормированные в 
Г,[с, 4]. Учитывая, что 4(5) Е Ну и 4\(5)Е Ну, при помощи 


разложения (6) получаем: 
м. й. 
р = ый = о 
п=1 


2 


р (8) 





№ 
Руа => (ое, а, + ди) 
п=1 


где ба", (5) = 4\, ($) —4у($). Отметим, что Су = 01" / о. 


Оценка для функционала И з Введем обозначения: 
04н(5) = 4(5)—-4(5); дих (х)=[Адах |(х); Ох =1-РБ - 


проектор из Ё.[а,Б] в ортогональное дополнение Н», пространства 
Ну. Из (7) следует, что |мл| > о |Чх|. Кроме того, Ом (5) Е Ну» и 


поэтому [би | < |АОх || 4х | . Воспользовавшись разложением 
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и(х) =их(х)-+ дин (х), из неравенства треугольника и предыдущих 
двух неравенств получаем оценку для У\? = |м Мин| | 
Ум <1+тй,, (М<М), (9) 
м = 2 
где ту = |АО»|/ о, Им = [64| 4х] = Им /\1 И» - 
Отметим, что оценка (9) завышена не более чем в м раз. В самом де- 
ле, существуют, очевидно, функци @4м(5), 04м(5), для которых 


[Мы = ом [Чм 


вольны, оценка 














бих | = АО || бах | Поскольку знаки этих функций произ- 


Й 


их + би | < || + [6 +2(их, би»), для них не мо- 
жет быть  априорно улучшена. — Менее точная — оценка 








[м ее би" < ([мх| Е [6и„ |)” для таких функций равносильна нера- 


венству (9); в то же время она отличается от предыдущей оценки не более 
чем в два раза, что без труда проверяется при произвольных функциях 


им (х), дим(Х). 


Оценки для функционала же Обозначим Оу = ХмЙм: Заметим, 
что 


п = ан - чу) + (а -@] 41 = а 14 +7. 
Отсюда и из (3) имеем [бам |< Хм бам | и, следовательно, 
[бах ак ры. Из (6) вытекает, что если 4(5)е Ну, то тогда 


м м 
оу || < |Ака|< 91 |4]. Используя эти соотношения, получаем две 
очевидные цепочки неравенств: 


о ПАчь а [Ау | -Аьб8 | _ 7 [88 











>1-Сурх! 
7’ |Ам9ы Ау»! дм |Ч» ыы 
гПа она ый 19] „ры. 
Ум” [Аух| ог || С 1Чм С» 


Таким образом, при не слишком больших значениях Рх справед- 
ливы оценки: 


< [у (Субх <1); (10) 
1- Сурм 

И о: | (р: (11) 
1 Ру 


Оценка (9) в сочетании с любой из оценок (10), (11) дает искомую 
оценку для функционала У» (5). При Ру >| функционал Я ‚ вообще 
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говоря, не является ограниченным сверху и может быть таковым лишь при 
неких специальных условиях. Подробный анализ случая Рим >] не пред- 


ставляет интереса, так как в приложениях величина Пу считается а рйой 


много меньшей единицы, а норма Жм должна быть меньше или порядка 


единицы (противное указывает на низкую эффективность проекционного 
метода и необходимость перейти к иной базисной системе [3]). 

Обсудим далее ряд вопросов, касающихся практического примене- 
ния полученных оценок (5), (9)-(11). 


Интерпретация параметра Ту и схема его численного расчета. По- 


строим пример, в котором значения Ту отыскиваются наиболее просто 
(все сопутствующие доказательства элементарны и для краткости опуска- 


НЕ, М ,. 
ются). Пусть {Ч (5)}„_! есть ортонормированный набор СФ оператора 


* Е = =: 
А А, отвечающих его СЧ д >09, >...б`и > 0. В каноническом методе 


проекций базисные СФ должны быть упорядочены по убыванию СЧ [1, 2]. 
тои версию проекционного алгоритма, в которой базисом 


У. (5) _, служат СФ, пронумерованные в ином порядке. При таком вы- 


боре базиа 7, =0, как и в канонической версии, а множество всех СЧ 


об, ‚ ге 1<п< М< М, совпадает с множеством {б`„ | — ‚ (среди СЧ 
С” встречаются одинаковые). В частности, может оказаться, что пере- 
становки не затрагивают СФ с номерами И < К. В этом случае ин 
Ту =бу, /Оу <1 при всех № < К. Если же в набор { }, (5)}. | вХО- 
дят СФ Ч, (5) с номерами т> М, то тогда Ту =б\, /б’,, где К - 
некие числа, такие, что К < п; при этом Ту, > 1. Отметим, что в данном 
примере У” =1, а точная оценка для функционала /‘\’ имеет вид 


1/2 
а < (1+ >) (ее максимальное отличие от оценки (9) в Я 


раз достигается при 777 = 1). В канонической версии метода проекций 


имеем Ту < |, что позволяет воспроизвести оценку из работы [6]: 
__ 2, 1/2 ти 


7м= = РИ) =а- пу) 
Е. и показывает, что в зависимости от свойств 
интегрального оператора и от выбранного базиса значения 7, могут быть 


как сколь угодно малыми, так и сколь угодно большими. При Ту в | 


1 
учет фактора УХ становится необходимым и существенно ограничивает 
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допустимые значения погрешности 77, при которых оценка (4) гаранти- 
и 
рует малую погрешность 77. 


Зависимость величин С, Ту от базиса может вызвать вопросы 
по поводу их минимизации. Заметим, однако, что при заданных значениях 
А, и и при отсутствии какой-либо иной априорной информации ника- 
кой метод решения интегрального уравнения (1) не может гарантировать 
реконструкцию проекции (5) с относительной погрешностью меньшей, 


1 

чем УЕ АЕ Последнее утверждение легко доказывается на основе раз- 

ложения (6). Оно отчасти очевидно, так как согласно (6) при любом выбо- 
м 

ре базиса СЧ О’, выступают весовыми коэффициентами, с которыми 


вклады от ортогональных проекций оригинала 9($) входят в правую 


часть интегрального уравнения. Как было показано, в отсутствие дополни- 
тельной информации оценка (9) не может быть значительно улучшена. 


Следовательно, значения С’,, Ту не характеризуют достоинства или не- 
достатки выбранной версии проекционного алгоритма, а отражают прин- 
ципиальные возможности реконструкции функций класса 4(5) Е Н ми 


функций, близких к элементам данного класса. (Единственной «мерой эф- 
фективности» проекционного метода, зависящей от базиса, является нор- 


ма Ух). 
В завершение настоящего раздела укажем удобный прием для чис- 
ленных расчетов параметра Ту. Воспользуемся — тождеством 
1/2 + 
[4о„|=|®»| ‚ де К, =(АОх) АОх. Вычисляя норму опера- 


тора №», можно считать, что он действует в пространстве Ну. Отсюда 


стандартным образом вытекает оценка 
1/4 


[40,|<| УВ |, (12) 


т,п=М-+1 
где В — (7-й) И Входящую В (12) сумму ряда можно оценивать 


численно, заменив в ней бесконечный верхний предел обрезающим пара- 
метром. 


Обозначим через ©'(М№, К’) наибольшее СЧ усеченной матрицы с 
элементами В, (М№М+1<т,п<К). Квадрат этого СЧ совпадает с 
нормой эрмитова оператора, получающегося из оператора №» при суже- 


нии его области определения до (К — №) - мерного подпространства, 
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ь к 
натянутого на систему функций {),(5)}„_х.!. Отсюда находим, что 


о (М, К) <| В, 








‚ или, что равносильно, 


(М, К) <|АОх |. (13) 
По принципу промежуточности о(М№, К+1) > о(М№,К), так что 


оценка (13) улучшается при увеличении А’. Последовательное увеличе- 
ние параметра К’ сопряжено с анализом плохо обусловленных матриц 
возрастающей размерности К-— №. Но наличие двусторонней оценки 
(12), (13) позволяет прервать этот процесс уже при сравнительно неболь- 
ших значениях А’. Число членов, удержанных в частичной сумме ряда из 
правой части неравенства (12), может варьироваться независимо от А’. 


Интерпретация функционала '/ ‘»). Множитель '/ > в оценке (4), (5) 
дает особого рода поправку связанную с — погрешностью 
04(5) =ах($) —4м(5). Эта погрешность возникает только при 


Ям = 0 и влияет на точность решений как непосредственно, так и кос- 


венным образом, увеличивая относительную невязку и тем самым понижая 
устойчивость решений к входным погрешностям. К сожалению, при 


С >>1 обе оценки (10), (11) сильно завышены, а получить в аналити- 


ческом виде более точные оценки не удается. 

При практической постановке некорректных обратных задач часто 
приходится анализировать перспективы реконструкции, не располагая до- 
статочной априорной информацией, в рамках неких качественных пред- 
ставлений, эвристических гипотез и т.п. Проекционный метод в основном 
как раз и предназначен для того, чтобы упростить и частично формализо- 
вать этот предварительный этап [4]. В связи с такими возможными прило- 
жениями обратим внимание на некоторые качественные свойства функци- 


онала и 

Заметим, что Оа(5) = Оубак (5$), где Ру:Ни— Ну - не- 
кий линейный оператор [3]. Допустим, что оператор Ру отображает про- 
странство Н № во все пространство Н ^ (это, как правило, выполняется, 
так кк АпаН м< Чт Н м= 00). Тогда ввиду произвольности функции 


4(5) и ее проекции Оу (5) функцию О (5) можно рассматривать как 


произвольный элемент пространства Ну. Функции 4м(5) и дах($) 
друг с другом никак не связаны, не считая ограничения 
[бам Мах| < ру. Далее предполагаем, что р, << 1. Из (7), (8) выте- 


2 
кает, что при таком условии большие значения У. >> 1 достигаются 





только при Су >> | и при одновременном выполнении следующих требо- 
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ваний: во-первых, функция 4/(5) должна быть квазиортогональна тем 
№ 

функциям Ч, (5), которые имеют малые номера И, т.е. отвечают боль- 

м 

шим СЧ О, ; 


доминировать в сумме (7); наконец, между функциями бах (5) и Чу (5) 


во-вторых, члены с такими номерами И должны, однако, 


(элементами векторного пространства Н ,) должна существовать опреде- 
ленная корреляция, при которой для указанных И выполнялись бы соот- 


м м 

ношения (Ф, ‚Ах + ба", ) << (Ф, ‚ан При введении некоторых 

гипотез о случайном выборе функций 4/(5), О (5) большие отклоне- 
2 

ния фактора и от единицы будут маловероятны. Это означает, что на 


2 
стадии эвристических прогнозов множителем Я в оценке (4) МОЖНО ПО- 


просту пренебречь. При непосредственном решении обратной задачи 
оценку (4) можно заменить апостериорными оценками [6], не содержащи- 


ми функционала Ух. 


Проиллюстрируем свойства функционала У” результатами чис- 
ленного эксперимента в духе метода Монте-Карло. Эксперимент был про- 
веден при значениях М№М=5; Я. =10 ро де п=12,..М; 
ру =0,1-—0,3. Функции ам(х), бам(х) моделировались двумя неза- 
висимыми случайными М - мерными векторами, ориентация которых отве- 


чает равномерному распределению в М- мерном кубе с центром в начале 
координат. Векторы нормировались из условия, что отношение 


|6 


тервале [0,Юх]. Определялись относительные частоты ©; наступления 


событий «у >», где #2 =1,2.3, й1=15; ИЕ, И5=2.5, 


В приведенной ниже таблице для каждой относительной частоты указан 
максимальный разброс ее значений, наблюдавшихся в 10 сериях из 10000 
испытаний. 





Ма» | есть случайная величина, равномерно распределенная в ин- 


Значения относительных частот ©; при различных значениях Рх 

















Юм о @> © 

0.1 0.0073 - 0.0103 0.0016 - 0.0026 0.0004 - 0.0015 
0.2 0.0329 - 0.0411 0.0102 - 0.0130 0.0045 - 0.0067 
0.3 0.0642 - 0.0717 0.0241 - 0.0277 0.0119 - 0.0162 
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Выводы. Недостающая оценка для функционала Ум получена путем его 


1 2 1 
факторизации на множители Уф ? ур . Анализ множителя Уф приводит 


к введению параметра 7, , который наряду с числом обусловленности 
определяет устойчивость решений к входным погрешностям. Множитель 


2 № 
и" отражает возможность резкого возрастания относительнои невязки 


интегрального уравнения при некой специфической корреляции между ор- 
тогональными проекциями функции-оригинала. При априорном исследова- 


о 2 № 
нии возможностей реконструкции фактор и" в порядке упрощающей 
гипотезы может игнорироваться. 
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\.М. ОВАСТЕЕУ 


ОМ ТНЕ ТМТЕСВАЕ ЕОЧАТТОМ ОЕТНЕ ЕТК$Т КТМО ВЕЗТОЦАЕ 
1М ТНЕ РКОЗЕСТТОМ МЕТНОО 


Гп {Пе ргемоцб рарегз Гог {Ле п\едга! РгедаНо!т едцаНоп оЁ {Пе Яг$Е Кта мии 
а дедепегае Кегпе! {Пе депегайтеа (поп гедиай2та) уег$юоп оЁ те ргодесЧоп 
те#о4 \м/аз Чеуеореа, Базе оп Ше агйгагу спосе оЁ Ле Баз!5 сотЫпеа 
МКР зоте а рйой езИтаНоп$ Гог {Пе еггог оЁ Пе зоивоп. ТПезе езбтаНоп$ 
дерепа оп Ме гезаца! мР/сП ма педесеЯ еа!ег. ТПе гези оР {Пе ргезепе 
рарег адтй {ю фаКе Ше гечаиа! по ассоипЕ апа $По\м/ Наб депегаЙу зреак- 
пд, {15 ассоипйпод 15 езепва|. 
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